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; Abstract 

> 

ly-j ■ The self- gravitating, spherically symmetric thin shells built of orbiting par- 

ticles are studied. Two new features are found. One is the minimal possible 
', value for an angular momentum of particles, above which elliptic orbits become 

possible. The second is the coexistence of both the wormhole solutions and the 
elliptic or hyperbolic orbits for the same values of the parameters (but differ- 
ent initial conditions). Possible applications of these results to astrophysics 
U" 1 and quantum black holes are briefly discussed. 



I. INTRODUCTION 



Thin shells play an important role in General Relativity. Because of the non-linearity of 
Einstein's equations we have to consider some simplified models that, describing essential 
features of the problem under consideration, can be solved analytically at least in quadrature. 
The thin shell formalism provides us with such a source of the gravitational field, the first 
integrals being incorporated in the very structure of the Einstein equations for the shells. 

This formalism was first developed in Jl| and applied to the gravitational collapse problem 
0. It was then almost forgotten for many years. Only in 1983 was it revived and used for 
studying dynamics of vacuum phase transitions in the Early Universe 0-10]. The possibility 



to obtain self-consistent solutions to the gravity equations together with the matter source 
equations has led to very interesting results (restrictions on the parameters of a vacuum 
decay process, enhancing of a vacuum decay probability around black holes and others) 
PJ6|J9]-|T2| , which were impossible to obtain by studying the dynamics of test particles (or 
fields) in a given background space-time. 

Comparing the self-consistent treatment of gravitating matter sources and the motion of 
the same source in a given gravitational field we should stress that the main advantage of 
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the former is that in this case we obtain an important information about a global space-time 
structure while in the latter case the matter fields "feel" space-time only locally. 

The very possibility to feel the space-time globally is extremely important in quantum 



theory. It was recently demonstrated |T^] in quantum black hole models by discovering a 
new quantum number. Thus, thin shell models serve as a powerful tool in the investigation 
of a variety of problems where effects of strong gravitational field may prove essential. 

In the present paper we studying a new class of spherically symmetric, self-gravitating 
thin shells. Namely, these shells are built of particles that follow elliptic or hyperbolic 
orbits, their pericentre being smeared uniformly upon the sphere with zero average angular 
momentum. The paper is organized as follows. In the Section II we briefly introduce the 
thin shell formalism with the emphasize on the spherically symmetric case. In the Section 
III the surface energy-momentum tensor is derived for the spherically symmetric shells with 
orbiting constituents. The Section IV is devoted to the dynamics of such shells. In the 
Section V we consider a special, very important case of light-like particles. And, finally, in 
the Section VI we discussed the obtained results and their application to astrophysics and 
quantum black holes. 

Throughout the paper we use the units % — c — 1 (h is the Planckian constant and c is 
the speed of light). In these units the Planckian mass rripi ~ 10~ 5 gr, length Ipi ~ 10~ 33 cm 
and time tpi ~ lCT 43 s are expressed in terms of the Newtonian constant G as G = nip 2 = 
1 pi — L pi- 



ll. GENERAL PRELIMINARIES 

1. Spherically symmetric space-times 

A general spherically symmetric space-time is a direct product of a two-dimensional 
sphere S 2 and a two-dimensional pseudo-Euclidean space-time M 2 . Therefore, a general 
line element can be written in the form 

ds = g^dy^dy" = r ) i kdx % dx k — R 2 {x)da 2 . (1) 

Here g^ v is a metric tensor of a four-dimensional space-time M 4 , y^ are the coordinates 
on M 4 , 7^ is a metric tensor of M 2 in a coordinate system x\ R(x) is a radius of a two- 
dimensional sphere and da 2 is the line element of S 2 , which can be parametrized with two 
angles, 9 and ip: 

da 2 = d6 2 + sin 2 9dip 2 . (2) 

The radius R(x) is defined in such a way that the area of the sphere equals 4irR 2 (x). 
Since Einstein's equations of General Relativity are invariant under general coordinate trans- 
formation we are able to subject altogether three components of two-dimensional metric 7^ 
to two conditions. Thus, a general spherically symmetric space-time can be locally described 
by only two functions: one is just the radius R(x) and another one comes from M 2 . It can 
be shown (for details see |9|,[T8|1) that for this second function we can choose an invariant 



Y k R,iRk, (3) 
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where j tk is the inverse to the metric tensor 7^, is a partial derivative of the radius with 
respect to a coordinate x 1 . The invariant function A is nothing but the square of the vector 
normal to the surfaces of constant radius R(x) =const. 

This normal vector can be space-like (A < 0), time-like (A > 0) or null (A = 0). It 
is simple to analyse these possibilities using orthogonal coordinates. Choosing some time 
coordinate t and some radial coordinate q, we have 

A = AR 2 — BR' 2 , A(t, q),B(t, q) > 0; (4) 

here dot and prime stand for derivatives along time (t) and radial (q) directions respectively. 
In the case of negative A the surfaces of constant radius are time-like. We have already used 
one of the two possible coordinate transformations (gauge freedom) to make coordinates 
orthogonal. Now we can exploit the remaining freedom to put q = R. The region where the 
radius can be chosen as a radial coordinate is called an i?-region. Similarly, for positive A 
and space-like surfaces of constant radius we can put t = R. The region where the radius 
can be chosen as a time coordinate is called a T-region. The notions of R- and T-regions 



were introduced in 15 



Further, in the .R-region (A < 0) we can never have R' = 0. So the sign of a partial 
derivative of the radius is an invariant. So, we have either R! > 0, which is called the 
i? + -region, or R! < 0, called the i?_-region. If we assume that a radial coordinate q ranges 
from —00 at the left to +00 at the right and call inner region the region to the left of the 
given surface R =const. and outer region the one to the right, then, in a _R + -region, the 
radii increase outside this given sphere, while they derease in the i?_-region. In the flat 
Minkowsky space-time A = —1, so we have globally an f? + -region and, can thus choose 
q = R everywhere. 

Analogously, in a T-region (A > 0) we can never have R = 0. Thus, there may be regions 
of inevitable expansions with R > 0, and regions of inevitable contraction with R < 0. The 
former are called T + -regions, while the latter are T_. 

The null surfaces of constant radius (A = 0) are called the apparent horizons and serve 
as boundaries between R- and T-regions. 

Thus, spherically symmetric space-time may have quite rich a structure. In general, it is 
some set of R± and T-t-regions separated by the apparent horizons. The most famous (and 
most important to us) example is the geodesically complete Schwarzschild space-time. The 
space-time is geodesically complete (or maximal analytically extended) if all the geodesies 
start and end either at infinities or at singularities. The Carter-Penrose conformal diagram 
for the Schwarzschild manifold is shown in Fig. 1 (on the conformal diagram all the infinities 
are brought to final distances). 

We see that in this simple case all possible regions are present. We have two isometric 
.R-regions (R_ to the left and R + to the right), each of them having an infinity, and two 
isometric T-regions (T + -region of inevitable expansion in the past and T_-region of inevitable 
contraction in the future, T + -region starting and T_-region ending with real singularities at 
R = 0, which are space-like surfaces). 

2. Thin shell formalism 

In General Relativity the inclusion of matter sources into considerations is not an easy 
task because of nonlinearity of Einstein's equation. Besides, the investigation of a self- 
consistent solution of a combined system of gravity and matter evolution is extremely im- 



3 



FIG. 1. Carter-Penrose diagram for the maximally extended Schwarzschild space-time. 



portant since it can differ substantially (and crucially) from the evolution of the same matter 
in a given background field. Thus there is a necessity in the construction of a very simple 
and tractable model for matter sources. One of these model are thin shells. The thin shells 
are viewed as a vanishing limit of thick shells of ordinary matter. In some cases, e.g. in 
cosmological phase transitions, such a limit can be justified 0. 

The general mathematical theory of thin shells in General Relativity was invented in [|TJ . 
It was formulated in a nice geometrical way as the set of equations for an extrinsic curvature 
tensor, which describes the embedding of a three-dimensional hyper-surface (world-sheet 
of a thin shell) into four-dimensional space-time (of course, a number of dimensions can 
be made arbitrary (d — 1) and d). Since a non-zero amount of energy is concentrated in 
a vanishing volume, such a thin shell is singular and, while all the metric coefficients are 
continuous across a singular hyper-surface, some of their derivatives undergo a jump. The 
equations for this jump are the following || 

[K{]-S{[Kl] = ^GSl (5) 

where K\ is the extrinsic curvature tensor describing the way of embedding of a time-like 
three-dimensional hyper-surface into a four-dimensional space-time (for a space-like hyper- 
surface one should change the sign of the right (or left) hand side). The square brackets 
denote a jump across a hyper-surface, [A] = A out — A in , Si is a surface energy momentum 
tensor of a shell. From the Einstein's equations just outside the shell there follow two more 
expressions 

Si\ j +[T?] = Q (6) 

{K}}S{ + [T:\ = 0. (7) 

Here {Kj} = |(-ftTj (out)— A'j(in)), and the index n denotes the outward direction normal to 
the shell. 
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In the spherically symmetric case everything becomes much easier. Of the total set of 
Eq.(B|) we are left with only two: 



[Ki] = AtcGSq (8) 

M + [Kl] = SnGSl (9) 
and from Eqs.(^|) and ([/]) we again have only two equations 

{K°}S° + 2{Kl)Sl + K] = (10) 

S° + ^(S o o -Si) + [TZ} = 0. (11) 

Here a dot means a derivative with respect to the proper time of the shell. Moreover, it 
appears |J, that Eq.(|l^) is an algebraic consequence of other equations, and Eq.(|ll|) can 
be considered as an integrability condition to Eqs.(||) and (^). Therefore, we can use any 
two of three equations, Eqs.(|]), @ and (|Tl|). Another advantage of spherical symmetry is 
that we are able to calculate a K|-component of the extrinsic curvature tensor for a general 
spherical space-time H. The result is 



Kl = -^?-A, (12) 

where p is the radius of the shell as a function of the proper time, p is its derivative, and A 
is an invariant function for a given spherical space-time, introduced earlier. What is cr? It 
is a sign function: a = +1 if the radii increase just outside the shell (in the outward normal 
direction) and a = +1 if they decrease. We know already that a = +1 in any i? + -region, 
and a = —1 in any R_ -region. The sign of a can change only in T_ -regions. 

At the end of this section we would like to make a rather important note. It concerns 
a qualitative difference between two approaches to the matter-gravitation dynamics. In 
the first approach one considers matter motion in a given background (say, Schwarzschild) 
metric. In the case of dust particles it means simply a geodesic motion. In the second case 
the motion is self-consistent, i.e. the full back-reaction of the matter fields on the space-time 
metric is taken into account. This is best illustrated by some simple model. Let us consider 
a spherically symmetric thin shell consisting of radially moving dust particles with total bare 
mass M; then Sg = Let them move in the gravitational field of a black hole of mass 

m. Owing to the back-reaction, the total Schwarzschild mass outside the shell will be, say, 
m + Am, Am > 0. Our main equation (M) is now 




2Gm , 2G(m + Am) GM 
aWp 2 + l ^ '-= , 13 

p V p p 

where a = a out (cr in = 1 because we put the shell outside the black hole). For the sake of 
simplicity we are interested only in a bound motion. Let po be a turning point, then 



2Gm / 2G(m + Am) GM 

aJl ^ J -= . (14) 

Po V Po Po 
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If cr = +1, the turning point is in the _R + -region of the outer metric and we call this a 
black hole case. For a = — 1, po lies in the i?_-region and we have a wormhole case. Let us 
introduce the following convenient notation e = jj and Ae = ^rp. It can be easily shown 
that 

a = +l if Ae > -(Ve^+I-e) (15) 



For e > 1 we have 



a = -l if Ae < -(Ve 2 + 1 - e). 



a = +l if Ae > — (16) 



a = -1 if Ae < 



1 

1 

4? 



In the Schwarzschild background limit Ae > £ we are left with only the i? + -region. Of 
course the shell can move in the Z?_-region as well, but in this limit the observer cannot 
distinguish between "in" and "out" (left and right on the Carter-Penrose diagram), so the 
shell can fill only local geometry. But with a self-consistent description we are able to probe 
a global geometry as well, and this may have far-going consequences for a quantisation 
procedure [[3 . 



III. THIN SHELLS WITH ORBITING CONSTITUENTS 

1. Construction. 

Our aim is to study the dynamics of a self-gravitating thin shells with orbiting con- 
stituents. What does the expression orbiting constituents mean? Let us consider a point 
mass moving according to law of Newtonian gravity (Kepler's problem in Celestial mechan- 
ics) in a field of gravitating centre. Its orbit is an ellipse with such a centre at one of the 
focuses. If we build an ensemble of such particles with the same angular momentum-to-mass 
ratio, they will have the same value of both pericentre and apocentre. Let us imagine that, 
initially, all these particles are smeared uniformly on a surface of a sphere whose radius 
is that of pericentre and let them start simultaneously (with equal absolute value of the 
velocity but in different directions, i.e. in different planes). Then such an ensemble will 
form a spherically symmetric thin shell oscillating between a pericentre and an apocentre. 
This is exactly what we call a thin shell with orbiting constituents. Such a construction can 
be applied equally to both non-relativistic and relativistic Coulomb problems. Of course, in 
relativistic case, orbits are no longer closed ellipses, but it is all the same qualitatively. 

We confine ourselves to the spherically symmetric case, not only because it is simple to 
treat. The equations for the gravitational field are non- linear, moreover, the back-reaction of 
the matter fields plays an important role in the case of strong gravitational fields, e.g. in the 
black hole physics we are mostly interested in, so that the addition of one source to another 
does not lead to a simple sum of the resulting gravitational fields. Besides, our distant goal is 
the construction of a quantum theory of such shells and we want to avoid the complications 
caused by the gravitational waves and, thus, encountering all the difficulties and obstacles of 
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the full quantum gravity. These unpleasant future are absent in the spherically symmetric 
models. 

The shells with orbiting constituents can also be useful in astrophysics, e.g. in studying 
dynamics of globular clusters. The very idea to construct spherically symmetric shells out 
of orbiting stars belongs to Bisnovatyi-Kogan But the full application of the idea 

to the analytical and numerical investigation of globular clusters requires consideration of 
intersecting shells, and this problem is far from being solve. 

2. Surface energy-momentum tensor. 

We already saw that, for a full description of the thin shell in the case of spherical sym- 
metry, we need only one equation, Eq.(|) (which is nothing but the energy conservation 
constraint). In order to solve it we should know the surface energy density Sq of our shell. 
This can be deduced from the integrability condition, Eq. (|TTD , provided we know an "equa- 
tion of state" Sf = S 2 (Sq) or that we calculate it independently. Before going to calculations, 
some general words are in order. 

The energy and the energy density depend, of course, on the choice of the time coordinate. 
In our case we need the surface energy density measured by an observer sitting on the shell 
(and not looking around). For this observer, the choice of the time coordinate is very poor, 
namely, t = t{r) where t is the proper time. But under such a reparametrization S$ is 
invariant (unlike 5*°° or Sqq). This invariance can also be understood in the following way. 
The angular components of the surface stress-energy tensor S 2 = S3 are invariant under 
transformations in a factorized two-dimensional space-time Mg. The trace Tr5? is also 
an invariant, and so is the mixed temporal component Sq —TrSj — S^ 2 ~ S% =TrSj — 2 >S 2 2 . 
Thus, Sq can depend only on the invariant functions of an underlying space-time (R and A). 
Moreover, Sq should not depend on A because in a two-dimensional space-time M 2 a world- 
sheet of our thin shell becomes a world-line, and by appropriate coordinate transformation 
any space-time can be made flat along the line (irrespective of the fact that Sq describes the 
internal structure of the shell and not its embedding properties). Thus, Sq depends only on 
the radius R of the shell and some appropriate integrals of motions. 

In order to get a surface energy density Sq we first calculate the sum of energies of the 
constituent particles in the frame of reference of an observer sitting on the shell. Then, 
taking into account that particles are uniformly distributed on the sphere of radius p we just 
divide the energy by Anp 2 . To fulfil this programme let us consider a point particle of mass 
/i. The action integral for such a particle is (see e.g. [I7j ) 

S = —fi J ds = (ds 2 = g a/3 dy a dy f3 ) = 

= -fxj \J g a p du a dvftds = -\l J \J g a p dy a dyP = (17) 

f / dy a dy 13 r , 

=~n r^"/ LdL 

Here u a = ^ is a four- velocity, t is some time coordinate (y° = t), and L is a Lagrangian 
and, thus, the energy heavily depends on the choice of time. Since we are interested in the 
spherically symmetric space-times seen by the observer sitting on the spherical shell the line 
element takes the form 

ds 2 = dr 2 -p 2 (r) (d9 2 + sin 2 6 dip 2 ), (18) 
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where p is the radius as a function of the proper time r. Then, as usual, we find the momenta 
p and the Hamiltonian (the dot means a derivative with respect to r): 



L = —p y 1 — p 2 2 — p 2 <p 2 , 



p p 2 sin 2 00 



pp 2 u e = —pue, 



'1 - p 2 9 2 - p 2 sin 2 OCp 2 

pp 2 sin 2 fly? 2 
£V = / = = PP u v = -pu v , (19) 

'l - p 2 2 -p 2 sin 2 Oip 2 



H = p e + p v tp - L = pu = 

= p Jl-u e u e -u^uv = \ P 2 + \pj + 2 1 2 q P. 
v V P p sin 6 1 
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Now, from the Hamiltonian equations 



• dH dH 

ope op^ 

dH 8H 

Pip 



10 dO 1 ^ dtp 
it follows that 

pI + -At^pI = const = j2 ( 21 ) 

sin r 

where J is the conserved angular momentum of each particle W71 . Hence, for the energy of 



our shell (remember that all particles have the same ratio a = -) we obtain 



E = M] jl + ^, (22) 

where M is the bare mass of the shell (i.e. the sum of all the rest mass of the constituent 
particles). Thus, for the surface energy density we have 



Having S® we can easily calculate S 2 = Sf. Indeed, the structure of the surface energy- 
momentum tensor S? copies that of the energy-momentum-tensor T£ for the dust matter 

Sf = Aui v? 

UiU^l, (24) 
except that now u l = since we are sitting on the shell. Then, we have 
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J2 fl 2 

u e u e + u ip u v = — = — -, 

[i A fr p z 

o « 2 
u u = 1 + — , 

P 2 



A 



M 



S 2 = S? 



Ma 2 



It can be readily checked that with these S® and S^ 2 the integrability condition inEq . (|TT|) is 
identically satisfied. 

Now we have everything for analysing a dynamical evolution of a self-gravitating thin 
shell with orbiting constituents. 

For further purposes we complete this section by writing the expression for Hamiltonian 
and energy for our thin shell in the specific background metric: 

ds 2 = Fdt 2 - -dr 2 - r 2 (d9 2 + sin 2 9 d^ 2 ) . (26) 
F 

This form includes the flat Minkovski space-time in Lorentzian coordinates for F = —A = 1, 
Schwarzschild space-time (i.e. a metric outside a spherically symmetric, neutral, massive 
object of mass m) for F = —A = 1 — 2hii» ) anc l Reissner-Nordstrom space-time (a metric 
outside a spherically symmetric, electrically charged, massive object of mass m and charge 
e) for F = —A = 1 — 2 ^ ni — — . The Hamiltonian is 

/ r r 



H = v / fJ/i 2 + Fp r 2 + ^. (27) 



Note the appearance of the radial momentum p r absent, for the on-shell observer. The 

dr' 
dt ■ 



expression for the shell energy written in terms of the radial velocity (^) reads as follows 



F 3 / 2 M a 2 
E= . =Wl + — . 28 



Note that now it is not the energy measured by an on-shell-observer, but the total energy 
(including all interactions a kinetic energies) measured by a distant observer (sitting near 
infinity). We can rewrite it in terms of the proper time: 

Fdt 2 _ I rfr 2 = dT 2 
j / , \ 2 / , \ 2 



F 



( dr 




[dt 


)'■(% 


dt\ 


2 1.2 


Tt) 


~F P - 
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(29) 




The last expression looks a little bit more elegant than Eq.(^). It is worth while to know 
that the same expressions for the energy follow from geodesic equations in the corresponding 
space-times, as it should be. 



IV. DYNAMICS OF THIN SHELLS WITH ORBITING PARTICLES 

In order to understand how General Relativity changes the dynamics of shells with 
orbiting particles, we first consider the well-known relativistic Coulomb problem, then we 
investigate the shells in the given Schwarzschild background and, at the last part of this 
Section, we show what is new that comes from accounting for a back-reaction of a shell 
dynamics on a space-time metric. 

1. Relativistic Coulomb problem 

To derive the main equation, which is just the energy conservation equation, we start 
with our general expression, Eq.(|M), for a charged thin shell with orbiting particles with 
Sq given by Eq . (p3|) . In this case the metric both inside and outside the shell is given, in 
general, by the Reissner-Nordstrom solution to the Einstein equation with 

F= _A=l-i^ + ^!, (30) 

where G is the Newton constant, m is the total mass (energy) of the system, and Q is its 
electric charge. We denote the mass and the charge inside the shell by m and Q, respectively, 
and that of the outside region by the (m + Am) and (Q + q), respectively. Then, the 
main equation reads as follows (in this subsection we need consider only the case with 

Pin — &out — +1) : 




2G(m + Am) G(Q + q) 2 



P 




(31) 



where M is the bare mass, a is the specific angular momentum of each particle, p(r) = r(t(r)) 
is the radius as a function of the proper time r on the shell, and dots denote the time 
derivatives. Taking now the limit of vanishing gravity (G — > 0) we obtain 




Am = dp* + 1. 1 + - + ^ ^ , (32) 



where Am is the total mass (energy) of the shell. For a hydrogen-like atom, Q = Z e, q = —e 
(electron charge), and after changing notation for the total energy (Am — > E) we have 
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v( P y 



a) 




FIG. 2. a) J < J cr , b) J > J a 



Ze 2 + \e 2 



(33) 



This last equation differs from the conventional form of the relativistic Coulomb problem in 
two respects. First, the kinetic term is different because, as we already explained, we use 
the proper time of the shell rather than the Minkowskian coordinate time. Then, there is 
an additional repulsive term proportional to e 2 , which describes the self- interaction of the 
charged spherical shell. Introducing a "potential" by Am(p = 0) = V(p) + 1, we get 



p 



(34) 



where by J cr — Z e — we denoted the critical value of the angular momentum. Remem- 
bering that a = -, we have in the case of a single orbiting particle (M = p) two different 
curves for our potential V(p) depending on whether J < J cr or J > J cr , as shown in Fig. 
2. 

2. Motion in the Schwarzschild background 

To investigate the shell motion in the given Schwarzschild background, we start with the 
same equation as in the preceding subsection, but with zero electric charge: 



2Gm 



(35) 



/„ 2Gm l 2G(m + Am) GM 
1 p 2 + 1 \ p 2 + 1 = i/l 



P 



P 



P 



P z 



The limit to a given background is the limit at which the total energy of the source E = Am 
is negligible with respect to the total energy m of the system. Therefore, we should expand 
Eq.(??) in Am up to the first order; the result is 



Am = E = MJp +1 - 



2Gm 



'1 + 



(36) 
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V(x) 


Ae-1 > 








2e 




Ae-1 < 




-1 









FIG. 3. 7 < 7cr = 2\/3e. The allowed regions are Ae > V(x) + 1. 



The formal difference from Eq.(|5BD for the relativistic Coulomb problem is that now the 
attractive force is incorporated into the kinetic term. Equation ( p6[) can also be considered 
as describing a shell consisting of radially moving particles with variable effective masses 
M e ff = M<Jl + s|. Introducing the dimensionless quantities Ae = ^p, e = jg, x = 
7 = and the "potential" V(x) = Ae(x = 0) — 1, we get 

Again, we have a critical value of the angular momentum 7 cr = 2a/3£ ( J cr = 2\/3G m /i, 
where \i is the mass of a single particle). The potential for 7 < 7 cr (J < J cr ) behaves 
qualitatively like the relativistic Coulomb potential in this case, apart from the fact that it 
now starts from x = 2 e (where it takes the value —1, the lowest possible value corresponding 
to Ae = 0). Thus, for small enough values of the angular momentum (7 < 7 cr ) we have two 
possibilities. Either Ae > 1, and the shell collapses starting from infinity (or, the other way 
around, it starts from the past singularity at p = and escapes to infinity). Or Ae > 1, 
and the shell starts from the past singularity at p = 0, reaches its maximal radius po at the 
turning point, and then recollapses to the future singularity at p = 0. 

When the angular momentum exceeds the critical value, the situation is essentially differ- 
ent from the case of the relativistic Coulomb problem. Starting from infinity, the potential 
first goes down to its minimum at x — |(1 — a/1 — £) e (£ = =-f-): 



VU = iy6 + |(l-(l-0 8/a )-l. (38) 
grows up to the maximum at x — 1(1 + \/l — £) e: 



V max = UQ + -Al + (l-0 3/2 )~h (39) 
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V(x) 


Ae-1 > 


2e 


-1 




\^<Ae-l<0 x 








/ " -1 < Ae-1 < V mtn 






FIG. 4. 2y/3e = j cr < 7 < 71 = 4e. 



and then again goes down to the minimal possible value (V = — 1) at the Schwarzschild 
horizon (xh — 2e). The origin of this new feature, the fall off of the potential near the 
Schwarzschild horizon is that in General Relativity any energy gravitates. So, in our case 
not only the bare mass M but also the energy of the angular momentum contribute to the 
gravitational attraction. This is reflected in the fact that the bare mass enters the equation 
only in the combination MJl + 

Figure 4 shows the behaviour of the potential in the case 2^/3e = 7 cr < 7 < 71 = 4e. 
Now we have already four possibilities. Case 1: Ae > 1. The unbound motion from infinity 
to the singularity (and vice versa). Case 2, V max + 1 < Ae < 1. This is the bound motion, 
the shell starts from the past singularity, expands up to the sufficiently large value of the 
radius and recollapses to the future singularity (analogously to the case Ae < 1 for 7 < 7 cr ) . 
In this case the shell has enough energy (Ae) to overcome the attraction due to a non-zero 
angular momentum. Case 3: V m i n + 1 < Ae < V max + 1. Here we have two possibilities 
depending on the initial conditions. First, we have orbits, as in the relativistic Coulomb 
problem. Second, we have bound motion with a collapse, but with much smaller value of 
the radius at the turning point (compared with case 2), and no analogue in the Coulomb 
problem. Case 4: < Ae < V min + 1. Only bound motion with a collapse. 

For 7 = 71 = 4e, the value of the potential at the maximum is just zero. Figure 5 
shows the potential in the case 7 > 71 = 4e. Again, we have four different cases. Case 
1: Ae > Vmax + 1 > 1. The total energy is large in order to overcome the repulsion due 
to the angular momentum, and we have a totally unbound motion (0 < x < 00). Case 2: 
1 < Ae < Vmax + 1- There are two possibilities, depending on the initial conditions. First, 
as in the Coulomb problem, the shell comes from infinity, reaches the turning point and 
goes back to infinity. Second, we have bound motion with a collapse, like for Ae < 1, (but 
now Ae > 1!). There is no such analogue in the Coulomb problem. And, last, Cases 3 and 
4 for Ae < 1 are qualitatively the same as Cases 3 and 4 for lower (7 cr < 7 < 71) angular 
momentum. 

At the end of this subsection we would like to note that from Eq.(^) it is easy to 
derive the radii for the smallest stable (7 = 7 cr ) and unstable (7 — ► 00, Ae = V max + 1) 
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FIG. 5. 7 > 7i = 4e. 

circular orbits. They are, of course, the same as for geodesies and equal to p s t = QGm and 
Punst = 3Gm, respectively. 
3. Self-gravitating shells 

Let us now turn to the most important case of self-gravitating shells with fill account 
for a back-reaction of the shell dynamics on the space-time structure. We start with the 
same equation as in the preceding subsection (all the notations are also the same) but with 
®out = & ( we will not consider here the shells Oi n = — 1 because there are no solutions in this 
case with the shell in the i? + -region which is most interesting for us). We get 




(40) 
and 

^f^Fi-hH)' (4i) 

while the potential reads now as follows: 

n^M-O)-.^^-!^)-!. (42, 

The new feature with respect to the case of a given background is the appearance of the 
self-interaction term in Eq. (f4ll) and, as a consequence, the possibility of changing the sign 
of a in Eq. (|40|) . This means that we can have not only the shells moving in the i? + -region 
(black hole case) but also those in the i?_-region (wormhole case). It is clear that in the 
wormhole case (a = —1) there cannot be unbound motion. The change of the sign of a and, 
hence, the transition from the black hole case to the wormhole case takes place when the 
turning point is exactly at the outer horizon, i.e. at Xq = 2 (e + Ae). It can be shown that, 
for fixed values of both ethere is exactly one value of 7 for which this condition holds. The 
other way around, for every e and 7 there is only one value of Ae a below which we have 
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a wormhole shell. This value of Ae a becomes negligibly small in the limit Ae <C e, i.e. in 
the Schwarzschild background limit, as it should be. So, we consider here another limiting 
case, e = 0, when inside the shell there are no other sources of gravitational field and the 
inner space-time is flat (in the general case some expressions can also be written but they 
are completely unreadable). Substituting for x the value x a = 2 (e + Ae) in Eq.(|40|) and 
putting x = we obtain 



Ae 1 



e + Ae x a \ 



1 + h (43) 



For e = it simplifies to 



1 

1 = —\ 



1 + ^, (44) 



with the following positive solution for Ae a : 



A£a = 271V 1 + v 1 + 4 ^ 2 - ( 45 ) 

If e 0, the corresponding value of Ae^ is less than for e = for given angular momentum 
7. The inequalities |^ > for the black hole case and |^ < for the wormhole case also 
hold. The new feature is the possibility for a wormhole shell to exist even for Ae > 1, which 
is completely forbidden for the shells with radially moving particles (there Ae < |)! 

Let us investigate now the behaviour of the potential V(x), Eq. (f4"2"D . The qualitative 
features are the same for all values of e, so that we will do this only for the simplest (and 
most tractable) case of zero inner mass, that is, we put e = 0. Thus, we have for the 
potential 



V = *e-l = Jl + £-±(l + £)-l. (46) 

Again, as in the case of the Schwarzschild background, there exists some critical value of the 
angular momentum (in our dimensionless notation it is 7 cr ) below which the potential looks 
as follows. It starts from the minimal possible value V = — 1 (corresponding to Ae = 0) 



at x m - m = T^TiY 1 + VI + I67 2 , crosses the horizontal line (Ae CT — 1) at i ff = 2Ae CT = 

+ \/l + 47" 2 ~, and then increases to zero at infinity. Below the line (Ae CT — 1) we have 
wormhole solutions, and above it there are only collapsing shells that start either from some 
turning point (for Ae CT < Ae < 1) or from infinity (for Ae > 1) at the _R + -region. 

To calculate the value of the critical angular momentum, we must solve (as in the 
Schwarzschild case) the following set of equations 

Ae 

V'(x) = = - - / = +TT^[1 + ^1 (47) 

r"u) ^_ == --[i + 6 Jl 

x d \ x z 
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FIG. 6. 7 < 7c 



Replacing the combination by £ (> 0) everywhere in the two last equations, we obtain 
a quadratic equation for £ with the positive solutions 



-3 + v/33 
12 ' 



(48) 



By a straightforward calculation we find 

(9 + v/33) 372 



8^ 



4.0859 



(49) 



2 _ (1 + V33)(9 + V33Y 

^Ycr 



3 • 2 7 



3.8184 



(50) 



Ae cr = 



(5 + v / 33)v / 9^v^3 = 



24 



0.8077. 



(51) 



For 7 > 7 cr the potential is qualitatively the same as in the Schwarzschild case, but now we 
may have coexisting wormhole shells and "ordinary" shells with the same parameters (but 
different initial conditions), the latter having either elliptic trajectories, as shown in Fig. 
7 or hyperbolic ones (Fig. 8) (note the line (Ae a — 1) never reaches the maximum of the 
potential). 

V. SHELLS WITH MASSLESS PARTICLES 



In this section we study the shells with orbiting massless particles. It should be stressed 
that althought particles are moving along null curves the shell itself is not null since a 
projection of a particle's velocity onto a radial coordinate line results in a time-like curve. 
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The corresponding limit is therefore very simple and straightforward. Indeed, let us look at 
our original equation (we confine ourselves to the a in = 1 case): 



/., , 2 Cm /.. , 2G(m + Am) GM I a 2 

— — — -=— 1 + 7 (52) 

and remember that the specific angular momentum a = ^ and M — N fj,, where J is an 
angular momentum of a single particle, fj, is its mass and N is a number of particles in the 
shell. Now, taking the limit \x — > we get 



. 9 , 2Gm I. 1 2G(m + Am) GL 
^ + — + l L_ U_ (53) 

Here we introduced a "total" angular momentum L = N //, although the real total angular 
momentum of our shell is, by construction, zero. Making the radius and masses dimension- 
less, p — Ipiy — VGy, m = m v [V = ^= (L is already dimensionless) we have 



Jy 2 + l aJy 2 + l-^ '- = — 54 

V y V y y 

Note that for y — ► oo the rapidity \y\ (a proper time velocity) tends to infinity as \y\ ~ 

\J ^fy, that is, the shell is accelerated to the speed of light. 

For the sake of simplicity, here we, again, consider only two limiting cases, namely, the 
limit of the Schwarzschild background Au <^ v, and the case v — (no other sources inside 
the shell). 

In the first case (Au u) we have 



Au = -Jy 2 + 1- — , (55) 

y V y 



and, introducing a potential U(y) = Au(y = 0), 



<W = if-T- (56) 

This potential has a maximum at y = 3 u with f/ max = ^r^;, as is shown in Fig. 9. We see 
that the higher the angular momentum L, the higher the maximal value U max , and the larger 
the black hole mass u, the lower the maximum. 

In the second case of inner mass u = we have for the turning points 



/ 2 Au L 

1 - a \ 1 - — = -2> ( 5? ) 



y V 2 y 2 / 

The changing in the sign of a occurs for y = 2 Au 



Au = ± (l - t%). (58) 
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The potential is shown in Fig. 10. Thus, above the maximum, the shells coming from 
infinity collapse; the less energetic shells have a reflexion point, i.e. they start to infinity, 
shrink up to the minimal radius and then accelerate back to infinity. At the same time there 
coexist shells with the same energy (but with different initial conditions) that start from the 
past singularities at p — 0, expand to the maximal radius, and then recollapse to the future 
singularity at p = 0. If their energy is above Av a , they live in the i? + -region (the black 
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hole case), while if the energy is lower they form wormholes. In the extreme limit Au = 0, 
such a wormhole closes up, and we have a closed Universe with only one thin shell with the 
constituent particles moving with the speed of light. 



VI. DISCUSSION 

In this paper we investigated the self-gravitating, spherically symmetric thin shells con- 
structed of orbiting particles. Is there anything new with respect to the well-known and 
well-understood geodesic motion in the Schwarzschild background? 

Yes, the self-consistent solution of Einstein's equations for a gravitational field, together 
with the equations for matter fields (dust particles in our case), leads to two new fatures. 

First, we found that there exists a minimal possible critical value for an angular momen- 
tum at which it becomes possible for particles to have an elliptic orbit. Such a minimum 
occurs for the shell inside which there is no other sources of gravitational field, and the inner 
space-time is, therefore, flat. In our dimensionless notation it is expressed as j cr ~ 1.95. In 
the conventional units it becomes 




J cr = N 7cr , (61) 



where mpi ~ 1CT 5 gr is the Planckian mass, p is the mass of an individual particle, and 
./V is the number of particles in the shell. Of course, this value is negligibly small for any 
atomic or subatomic system. But for the other extreme situation, say, for stellar nebula, it 
may be rather important. Indeed, the stars in this case moving along elliptic orbits around 
a common centre of mass can be approximated by a set of "thin" spherically symmetric 
shells with orbiting dust particles. For the innermost layer we have just the minimal critical 
angular momentum. The ratio (^j-) is now huge, as the number of stars could be. But how 
about the value J cr ? Is it large or small? To understand this let us calculate the angular 
velocity uj of the particles moving along circular orbit with J = J cr and radius p — p, 
Using the expression for the four-velocities at the very end of Section III, we easily obtain 



cr ■ 



Pi y/l + > 

a cr = — . (62) 
/i 

Thus, the linear velocity V does not depend on both mass p and number N of particles and 
reads 

V- = a r „ = ^— = I 1 7=r 1 - o I. (63) 



Per yjl + ^ X crX fl + 



-y 2 

I cr 




It can be shown that the linear velocity of the circular orbits (at the minimum of the 
potential) goes down with increasing 7. So, the maximal possible linear velocity for the 
stable circular orbit is 40% of the speed of light (the latter is just 1 in the chosen units). 
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This limit should be compared with the corresponding limit for a geodesic motion in the 
Schwarzschild background, where — = 4? and v CT = \. Thus, the account for back-reaction 
of the matter fields on the gravitational field leads in our case to a more severe restriction of 
the maximal possible velocity for the circular orbits. Maybe this will help us to understand 
how black holes could be formed in the galactic nuclei. 

The second new features is the wormhole solutions. The necessity to take into account 
these (rather unusual) solutions becomes crucial in quantum theory. It was shown in |14| that 
the very existence of the wormhole R- -region requires a new quantum number for "good" 
wave functions. For example, for simplest bound motions, we need two quantum numbers 
(with only one in the conventional quantum mechanics), and, what is more surprising, we 
have a quantum number for an unbound motion, in particular for radiating light-like (null) 
shells. It seems that it is a radiation spectrum that is responsible for a quantum black hole 
mass spectrum; and structure of matter field energy levels is responsible for the black hole 
entropy. To get a highly degenerate black hole mass spectrum (and, thus, a huge entropy) we 
need quite rich a structure of both radiation and matter source spectra. The features of the 
potential considered in this paper seem sufficiently rich for such purposes. We mean that in 
the case of the orbiting particles we have not only the existence of the wormhole regions but 
also the coexistence of both the elliptic or hyperbolic orbits and wormhole solutions for the 
same values of the parameters (but different initial conditions). Such a coexistence will lead 
to a double-splitting of the energy levels as in the famous quantum mechanical double-well 
problem. 
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